Spin-Charge Coupling in lightly doped Nd 2 - x Ce x Cu04 
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We use neutron scattering to study the influence of a magnetic field on spin structures of Nd2Cu04. 
On cooling from room temperature, Nd2Cu04 goes through a series of antiferromagnetic (AF) phase 
transitions with different noncollinear spin structures. While a c-axis aligned magnetic field does 
not alter the basic zero-field noncollinear spin structures, a field parallel to the Cu02 plane can 
transform the noncollinear structure to a collinear one ("spin-flop" transition), induce magnetic 
disorder along the c-axis, and cause hysteresis in the AF phase transitions. By comparing these 
results directly to the magnetoresistance (MR) measurements of Ndi.gTsCeo.ccsCuCU, which has 
essentially the same AF structures as Nd2Cu04, we find that a magnetic-field-induced spin-flop 
transition, AF phase hysteresis, and spin c-axis disorder all affect the transport properties of the 
material. Our results thus provide direct evidence for the existence of a strong spin-charge coupling 
in electron-doped copper oxides. 

PACS numbers: 74.25.Fy, 74.72. Jt, 75.25.+Z 



I. INTRODUCTION 

Understanding the role of magnetism in the trans- 
port properties and superconductivity of high-transition- 
temperature (high-T c ) copper oxides remains one of the 
important unresolved problems in the physics of transi- 
tion metal oxidesi. The parent compounds of high-T c 
cuprates are antiferromagnetic (AF) ordered Mott insu- 
lators composed of two-dimensional (2D) CuC>2 planes. 
When holes or electrons are doped into these planes, the 
long-range AF ordered phase is destroyed and the cop- 
per oxide materials become metallic and superconduct- 
ing with persistent short-range AF spin correlations 2 *^. 
While much work over the past decade has focused on 
the interplay between magnetism and superconductivity 
because spin fluctuations may mediate electron pairing 
for superconductivity^, understanding the relationship 
between AF order and transport properties through the 
metal-insulator transitions (MIT) in these doped copper 
oxides is interesting in its own right. 

For instance, the parent compounds of hole-doped 



cuprates have collinear AF spin structure, where each 
Cu 2+ spin is aligned opposite to its neighbors^. For 
La2Cu04 and La2- a; Sr a ;Cu04 in the lightly doped re- 
gion, the Cu 2+ spins in the CuC>2 planes are slightly 
canted from the direction of the staggered magnetiza- 
tion to form a weak ferromagnetic (FM) momenljL&ii. 
As a consequence, an applied external magnetic field 
can manipulate the AF domain structure and induce 
large anisotropic magnetoresistance (MR) effect 10 . In the 
case of the parent compounds of electron-doped materials 
such as Nd2Cu04 and Pr2Cu04, the magnetic structures 
are noncollinear, where spins in adjacent CuC>2 layers are 
90 degrees (90°) from each other [Fig. 1(c)], due to the 
pseudo-dipolar interaction between the rare-earth (Nd 3+ 
and Pr 3+ ) and Cu 2+ ionsiiiiSiH. Application of a mag- 
netic field in the CuC>2 planes will induce a "spin-flop" 
transition by transforming the noncollinear structure to a 
collinear onei^i^i^, and the critical field (Bsf) depends 
on the direction of the magnetic field with the field along 
the Cu-Cu (B| | [110]) direction generally has a smaller 
Bgi^ii. A c-axis aligned magnetic field has no effect on 
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the noncollinear spin structure— 

Recently, Lavrov et al*& have reported that 
an in-plane magnetic field can induce a large 
MR effect in lightly electron-doped copper oxide 
Pri. 3 _ :r Lao.7Ce x Cu04 (x — 0.01). The authors find 
four-fold-symmetric angular-dependent MR oscillations 
for Pri.2gLao.7Ceo.oiCu04 in the low-temperature non- 
metallic regime. Similar data have been also obtained in 
nonsuperconducting Pri.gsCeo.isCuCU independently— 
Since Pri.2gLao.7Ceo.oiCu04 has a noncollinear spin 
structure at low temperatures and the critical fields for 
spin-flop transition and MR effects are similar, the novel 
MR phenomenon in this material has been attributed to 
the spin structure rearrangement from the noncollinear to 
the collinear stated. For a c-axis aligned magnetic field, 
the observed negative MR effect in the normal state of 
several different electron-doped cuprates has been inter- 
preted as a result of two dimensional weak localization 
by disorder—, Kondo scattering from Cu 2+ spins in the 
CuC>2 planed, or spin scattering from field-induced mag- 
netic droplets formed around impurities— 

While these recent MR measurements on 
Pri.2gLao.7Ceo.oiCu04 and Pri.ssCeo.isCuCU clearly 
suggest a close coupling between spin-flop transition 
and MR effects^, the data may also be interpreted 
as partial rearrangement of magnetic domain walls by 
magnetic field to allow conductivity of electrons along 
a preferred direction—. In the latter case, the magnetic 
domain walls are segregated by the doped charge 
carriers into inhomogeneous patterns, such as stripes— 
If the transport properties in lightly electron-doped 
cuprates are indeed determined by spin reorientations 
and not by stripes, one would expect intimate correla- 
tions between the MR effects and spin-flop transitions 
in other families of electron-doped materials. Since 
Nd2CuC>4 exhibits three AF phase transitions with dif- 
ferent noncollinear spin structures on cooling from room 
temperaturei^i^*i^, far more complicated than the single 
AF phase transition found in Pri^gLao^Ceo.oiCuO^i 9 - or 
Pri.ssCeo.isCuO^Sil, a combined neutron scattering and 
MR investigation should shed new light on the interplay 
between spin and charge coupling in the material. 

In this article, we describe our neutron scattering and 
MR measurements on single crystals of Nd2CuC>4 and 
lightly electron-doped Ndi.g75Ceo.o2sCu04, respectively. 
For neutron scattering, we choose to study Nd2CuC>4 be- 
cause of its complicated AF phase transitions 14 . While 
previous work showed that a magnetic field applied par- 
allel to the CuC>2 planes transforms the spins from 
the noncollinear to collinear AF structur e - 14 ' 15 , there is 
no systematic work on how the field-induced collinear 
spin structure affects the zero-field AF phase transi- 
tions. We find that application of a B||[110] field can 
induce c-axis spin disorder and hysteresis in the AF 
phase transitions. Since lightly electron-doping the in- 
sulating Nd2CuC>4 induces enough charge carriers to al- 
low transport measurements 19 but does not change its 
basic AF spin structures^, we compare the MR effects 
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FIG. 1: NcbCuCU spin structures in (a) type-I/III and (b) 
type-II noncollinear states, where spins are indicated by the 
arrows. X, Y and Z represent the interactions between 
Nd-Nd, Nd-Cu and Cu-Cu spins respectively, as defined by 
Sachidanandam et alw—. (c) the schematic phase diagram of 
CUO2 planes in different phases at zero-field (bottom row) 
and Bsf||[110] (top row). Here only Cu spins are shown for 
clarity. The filled and unfilled circles represent L — and 1 /2 
layers of Cu atoms respectively. The Ti , T2 , and T3 represent 
transition temperatures for the three different noncollinear 
phases at zero fieldii. 



in Ndi.g75Ceo.o25Cu04 to the neutron scattering results 
on Nd2CuC>4. Surprisingly, we find that the transport 
properties of Ndi.g75Ceo.o2sCu04 are very sensitive to 
the modifications of spin structures in the system. Our 
results thus provide further evidence for the existence of 
a strong spin-charge coupling in electron-doped copper 
oxides. The organization of this article is as follows. In 
Section II, we describe the experimental setup for neu- 
tron scattering and transport measurements. Our neu- 
tron scattering results are presented in Section III while 
MR transport data are shown in Section IV. In Section V, 
we compare the neutron scattering and transport data. 
Finally, Section VI summarizes the conclusions of our 
work. 
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II. EXPERIMENTAL SETUP 



We grew single crystals of NCI2CUO4 and 
Ndi.97 5 Ce .o25CuC>4 by the traveling solvent floating- 
zone method. The samples were grown at a speed of 1 
mm/hour under 4 atm O2 pressure in a sealed quartz 
tube 26 . All the crystals are single domain as confirmed 
by a polarizing light microscope and Laue X-ray diffrac- 
tion. The Nd2CuC>4 single crystals are cylindrical and 
have dimensions of about 4 mm in diameter and 15 mm 
in length. 

The neutron scattering measurements on Nd2CuC>4 
were performed on the HB-1 and HB-3 triple-axis 
spectrometers at the high- flux-isotope reactor (HFIR), 
Oak Ridge National Laboratory (ORNL). We specify 
the momentum transfer (q x ,qy,qz) in units of A -1 as 
(H,K,L) = (q x a / 2-rt , q y b / 2tt , q z c / 2-rt) in reciprocal lat- 
tice units (r.l.u.). The lattice parameters of the tetrag- 
onal unit cells of Nd2Cu04 are a = b = 3.944 A and 
c = 12.169 A. To prevent the samples from rotating un- 
der the influence of a magnetic field, they were clamped 
on solid aluminum brackets and placed inside a 7-T verti- 
cal field magneti 1 ^. For the experiment, we use pyrolytic 
graphite as monochromator, analyzer and filters. The 
collimations were, proceeding from the reactor to the de- 
tector, 48'-40'-sample-40'-120' (full width at half maxi- 
mum or FWHM), and the final neutron energy was fixed 
at Ef = 14.78 meV. The experiments were performed in 
the [H, H, L] scattering plane where the applied vertical 
field is along the [110] direction (B||[110]) in the Cu0 2 
plane. 

For our transport studies, we align and cut one 
Ndi.975Ceo.o25Cu04 single crystal into three rectangu- 
lar blocks, whose typical size was 3x1x0.5 mm 3 . Using 
the regular four-points method, as shown in the inset of 
Fig. 2, the resistance was measured by the AC-transport 
option of a commercial 14-T physical property measure- 
ment system (PPMS). The AC currents were chosen 
along the [100], [110], and [001] directions and the corre- 
sponding resistance were labeled as -Rnoo] , -R[no] > an d R c , 
respectively. Since Ndi.97sCeo.o25Cu04 has tetragonal 
crystal structure, the a- and 6-axes are indistinguishable. 
As a consequence, a magnetic field along the [100] ([110]) 
direction is equivalent to the [010] ([110]/[110]) direction. 
Figure 2 shows the temperature dependence of zero-field 
resistances in the three high symmetry directions. Simi- 
lar to lightly doped Pri.29Lao.7Ceo.oiCu04^, resistances 
in all directions of Ndi.g75Ceo.o25Cu04 increase with de- 
creasing temperature and show an insulating behavior 
at low temperatures. In addition, the resistivity data 
show no indication of the influence of AF phase transi- 
tions. Since Nd1.975Ceo.025Cu.O4 is only slightly doped 
away from Nd2Cu04, it is reasonable to assume that 
these two systems have similar spin structures and AF 
phase transitions^. 
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FIG. 2: The temperature dependence of R c , i2[ioo] and R[no] 
at zero field, where the subscripts represent the current di- 
rections. The resistances in all three directions go up rapidly 
with decreasing temperature. The inset shows the regular 
four-points setup. Throughout the measurement, the c-axis 
of the crystal is always along the axis of rotation and the 
magnetic field rotates within the CUO2 (ab) plane. 

III. NEUTRON SCATTERING RESULTS 

Before describing our results on the influence of an in- 
plane magnetic field on the magnetic order of Nd2Cu04, 
we briefly review its zero-field spin structure a 14 ' 15 . When 
Nd2Cu04 is cooled from room temperature to T\ « 275 
K, its Cu spins first order into the noncollinear type-I 
spin structure of Figs. 1 (a) and 1(c). On further cooling 
to T 2 ~ 75 K, the Cu spins in the adjacent layer rotate by 
180° about the c-axis from the type-I phase and reorient 
the system into type-II phase [Fig. 1(b)]. Finally below 
T3 k 30 K, the Cu spins rotate back to their original di- 
rection to form the type-Ill phase [see Fig. l(c)}i£. Mag- 
netic structure factors (|F(l/2, 1/2, L)\) for noncollinear 
type-I and III phases at (1/2, 1/2, L) positions are 

\F(L = odd) | 2 = 32( 7 e 2 /2mc 2 ) 2 

\fcuMcu + 2cos(2irLz)f Nd M Nd \ 2 ; 

\F(L = even) | 2 = 32( 7 e 2 /2mc 2 ) 2 

(2aL/c) 2 /(2 + (2aL/c) 2 ) 
\fcuMcu + 2cos{2irLz)f Nd M Nd \ 2 . 

(1) 

For type-II noncollinear spin structure, we have 

\F(L = odd) | 2 = 32( 7 e 2 /2mc 2 ) 2 

(2aL/c) 2 /(2 + (2aL/c) 2 ) 
\fcuMcu + 2cos(27rLz)f Nd M Nd \ 2 ; 
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\F(L 



2\2 



32(7e 2 /2mc 2 ) 

\fcuMcu 



2 cos(2tt Lz)f Nd M Nd \ 2 ; 



(2) 



where je 2 /2mc 2 = 0.2695 x 1(T 12 cm, z = 0.35, f Cu , 
fNd, Mcu, and Mcu are magnetic form factors and or- 
dered magnetic moments for Cu and Nd ions, respec- 
tively. 

From the magnetic structure factor calculations, we 
find that the intensities of AF Bragg peaks at the 
(1/2, 1/2, L) positions depend sensitively on the detailed 
spins arrangement. For example, when Nd2CuC>4 is 
cooled below the spins reorientation transition temper- 
ature T 2 [Fig. 1(c)], intensity of (1/2,1/2,1) peak de- 
creases while that of (1/2, 1/2, 2) increases^. On further 
cooling to below T 3 , the scattering at these positions re- 
cover to their high temperature values [Figs. 3(a) and 
3(b)]. Gaussian fits to the data show that the scatter- 
ing is resolution limited with the FWHM of 0.02 r.l.u 
for (1/2,1/2,1) (or AL = 0.02 r.l.u.). In principle, one 
should calculate the coherence length of a Bragg peak 
from the formula of iV-slit grating diffraction^. How- 
ever, the lincshape of our observed diffraction peaks is 
well described by a Gaussian, equivalent to the iV-slit 
function in the limit of large N. By Fourier transform of 
the Gaussian peak, we estimate a minimum spin-spin co- 
herence length of -530 A using CL = [41n(2)/7r](c/AL). 

When a 5-T magnetic field is applied along the [1,1,0] 
direction, the noncollinear spin structures at different 
temperatures are transformed into collinear spin struc- 
tures [Fig. 1(c)]. In phases I and III, (1/2, 1/2, L = odd) 
peaks vanish while (1/2, 1/2, L = even) are enhanced 
with magnetic structure factors as 



\Fc(L 



64( 7 e 2 /2mc 2 ) 2 
(2aL/c) 2 /(2+ {2aL/c) 2 ) 
\fcuM Cu + 2 cos(2n Lz)f Nd M Nd \ 



(3) 

For phase II, (1/2, 1/2, L = even) peaks vanish while 
(1/2, 1/2, L = odd) reflections change as 

\F C (L = odd)| 2 = 64( 7 e 2 /2mc 2 ) 2 

\fcuMcu + 2cos{2irLz)f Nd M Nd \ 2 . 

(4) 

Figures 3(c) and 3(d) confirm that (1/2,1/2,1) is en- 
hanced while (1/2, 1/2, 2) vanishes after spin- flop transi- 
tion. Figure 3(e) shows the integrated intensities at the 
(1/2,1/2,1) and (1/2,1/2,2) positions as a function of 
increasing temperature after a 5-T magnetic field is ap- 
plied at 15 K to induce the type-Ill collinear state [Fig. 
1(c)]. From the temperature dependence of their inten- 
sities, it is clear that spin-flop transitions occur between 
the three collinear states at similar temperatures as that 
of the zero-field AF phase transitions^. 

To test whether the field-induced type-II collinear spin 
structure depends on the magnetic field hysteresis, we 
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FIG. 3: Typical neutron scattering results around 
(1/2, 1/2, 1) and (1/2, 1/2,2) Bragg peaks in the (a)(b) B = 
spin noncollinear and (c)(d) B = 5 Tesla spin collinear states. 
We probed the magnetic Bragg peaks at 15 K, 40 K, and 80 K 
in type III, II, and I phases, respectively. All peaks are fit by 
Gaussians on sloped backgrounds as shown by the solid lines, 
(e) The temperature dependence of the integrated intensities 
at (1/2, 1/2, 1) and (1/2, 1/2, 2) positions at B = 5 Tesla. The 
type-Ill to II and type-II to I collinear phase transitions are 
around 30 K and 70 K, respectively. 



perform neutron experiments at 40 K in two ways. We 
first cool the sample at zero-field to 40 K and then in- 
crease the B| | [110] field to 5-T as sketched in the inset 
of Fig. 4(a) (process a). The [1/2, 1/2, L] scan from 
0.8 < L < 5.2 r.l.u. shows resolution-limited peaks 
around L = 1,3 and 5 as shown in Fig. 4(b). Since the 
(1/2, 1/2, 1) peak has a Gaussian lincshape with FWHM 
of 0.02 r.l.u. along the L direction, we estimate that the 
c-axis magnetic coherence length is around 530 A. Now, 
if we applied the 5-T field at 15 K and then increased the 
temperature to 40 K (process b), the (1/2, 1/2, L) (L = 
odd) peaks remain but with much different lineshape 
[Fig. 4(c)]. Instead of having resolution-limited Gaus- 
sian lineshape, the peaks are Lorentzians (The Ornstein- 
Zernike form) with much broader widths but the same 
integrated intensity as that in Fig. 4(b). For example, 
while the FWHM of the (1/2,1/2,1) peak increases to 
0.046 r.l.u. from 0.02, its peak intensity also drops by 
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FIG. 4: (a) The in-plane [H,H,1] scans across (1/2,1/2,1) 
in two different processes at 40 K. The schematic diagrams 
of processes a and b are shown in the inset, (b) The long 
[1/2, 1/2, L] scan in process a, where L changes from 0.8 to 5.2 
r.l.u. Only three resolution-limited peaks are found around 
L = 1, 3 and 5. (c) Same scans as (b), but through process b. 
Note that the L-widths of these peaks become considerably 
broader. The temperature and field in (b) are identical as 
that of (c), i.e. T — 40 K and B = 5 Tesla. However, the 
processes of applying field are different, as shown in the inset 
of (a). The peaks in (b) and (c) are fitted by Gaussians and 
Lorentzians, respectively. 



^50% [Figs. 4(b) and 4(c)]. These observations suggest 
that the c-axis spin-spin correlation function decays ex- 
ponentially with a much shorter coherence length [^84 A 
using 1/k where k is HWHM in A -1 of the Lorentzian in 
Fig. 4(c)] in process b. Since the only difference between 
the type-Ill and II collinear states is the 180° spins ro- 
tation in adjacent layers [Fig. 1(c)], the short spin-spin 
coherence length in Fig. 4(c) suggest the presence of a 
c-axis spin disorder. On the other hand, in-plane scans 
along the [H, H, L] (L = 1, 3, 5) direction only show 
slight broadening when field is applied at low tempera- 
ture in process b [Fig. 4(a)], thus suggesting most of the 
spin disorder occurs along the c-axis. 

To further investigate how hysteresis in application 
of a magnetic field can affect the spin arrangements in 
Nd 2 Cu04, we studied its c-axis coherence lengths in all 
three collinear phases shown in Fig. 1(c). There are three 
different ways to reach the expected temperature and 
field of 40 K and 5-T in the type-II collinear spin phase, 
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FIG. 5: Magnetic field hysteresis effects on AF spin structures 
of Nd2Cu04 in all three phases, (a) (1/2, 1/2, 1) Bragg po- 
sition in phase II and (b),(c) (1/2, 1/2,2) position in phase I 
and III, respetively. The corresponding processes of field and 
temperature are shown in the inset of each figure. 



which arc labeled as a, b and c in the inset of Fig. 5(a). 
Precesses a and b are the same as Fig. 4 and process c 
involves getting to 80 K in zero field, applying the 5-T 
field, and then cooling the sample to 40 K. In other words, 
the spin system changes from the noncollinear state to 
the collinear state and remains in type-II collinear spin 
phase during the process a, whereas the system under- 
goes phase transitions from type-Ill and type-I collinear 
states to type-II collinear state in processes b and c, re- 
spectively [Fig. 1(c)]. 

By comparing the (1/2, 1/2, 1) data of Fig. 4 in pro- 
cesses a and b with that of c in Fig. 5(a), it becomes 
clear that processes b and c have considerably broader 
widths. This means that b and c processes induce large 
c-axis spins disorder and the amount of disorder depends 
on the history of field application. When the system is in 
type-I phase at 80 K, the collinear spin phase can be in- 
duced either by simply applying a field at 80 K [process c 
in the inset of Fig. 5(b)], or by having a type-II collinear 
phase at 40 K and then increasing the temperature to 
80 K [process a in the inset of Fig. 5(b)]. Clearly, the 
[1/2, 1/2, L] (L — 2) scans in Fig. 5(b) show resolution- 
limited Gaussian peak in c but two components lineshape 
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with a sharp Gaussian peak on a broad Lorentzian back- 
ground in process a. These results again suggest that 
magnetic field hysteresis affects mostly the spin arrange- 
ments along the c-axis. Finally, when temperature is in 
the type-Ill collinear phase at 15 K, we find no obvious 
difference between the two processes in the inset of Fig. 
5(c): one of which is applying field at 40 K and then 
decreasing the temperature to 15 K (process a), and the 
other is applying the field at 15 K directly (process b). 
Therefore, spin disorder appears whenever phase transi- 
tion occurs between two different collinear states except 
for the transition from type-II to III. 

Figure 6 shows the magnetic field dependence of the 
(1/2,1/2,1) peak at T = 40 K under different condi- 
tions. When a B||[110] field is applied from the zero-field 
noncollinear type-II state at 40 K [process a in the in- 
set of Fig. 6(a)], a noncollinear to collinear spin-flop 
transition occurs around 1-T, and the FWHM of the 
peak does not change during the process [Figs. 6(a) 
and 6(b)]. This suggests that the entire AF structure 
responds to the influence of the applied field. If we warm 
to 40 K using process b shown in the inset of Fig. 6(a), 
the (1/2, 1/2, 1) peak has a broad FWHM along the c- 
axis but with the same integrated intensity as process a. 
As a function of decreasing magnetic field at 40 K, the 
FWHM of (1/2, 1/2, L) at L = 1 decreases continuously 
until reaching the value of process a. This suggests that 
long-range c-axis spin coherence length is restored in the 
process. 

Because the spin disorder is related to the phase tran- 
sitions under field, it is natural to ask what happens to 
the transition itself under different conditions. Here, we 
use two ways to pass the type-II to type-Ill phase tran- 
sition, and label them as processes a and d in the inset 
of Fig. 7(a). We carefully monitor the (1/2, 1/2, 1) peak 
as the system passing the transition under different con- 
ditions. At each measured temperature, we wait 10 to 
15 minutes to ensure that the peak intensity has no time 
dependence and the outcome of the experiment is sum- 
marized in Fig. 7. Surprisingly, there is a clear hysteresis 
in the phase transition behavior, i.e. the transition tem- 
perature of process b is about 2 K higher than that of 
process a. 

To understand how this hysteresis occurs, we decreased 
the temperature following process b, i.e. we probe the 
(1/2, 1/2, 1) peak in process d after b. Figure 7(a) indi- 
cates that the transition temperatures of processes a and 
d are identical. While this suggests that spin disorder 
itself may not induce the lower transition temperature in 
decreasing temperature processes, the observed hystere- 
sis must be associated with the field-induced collinear 
spin structures and their free energy differences in dif- 
ferent phases. Similar hysteresis behavior can also be 
found for the transition between collinear type-II and I 
phases [Fig. 7(b)]. In this case, the integrated intensities 
of (1/2,1/2,1) during process c are much smaller than 
those during process a below 70 K, while the intensi- 
ties of (1/2, 1/2, 2) remain zero thus ruling out a possible 
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mixture of the two phases. 

In previous work on Pr2Cu04 2§ , which has a non- 
collinear type-I/III spin structure [Fig. l(a)Ji£, diffuse 
scattering associated with interplane short-range order 
was observed above the spin-flop transition critical field 
(3. 1-T) around the forbidden AF position (1/2,1/2,1) 
at T = 1.5 K. Since similar diffuse scattering was not 
observed around (1/2,1/2,2), the authors suggest that 
the diffuse scattering arises from the persistent mid- 
range interplane correlations^. For Nd2Cu04, we find 
no evidence for similar short-range order spin order at 
(1/2,1/2,1) (Figs. 2-5). Instead, we find clear evidence 
for field-induced AF phase transition hysteresis and spin 
disorder. If the spin degree of freedom is strongly cou- 
pled to the charge carriers, one would expect to observe 
changes in electrical transport properties uniquely asso- 
ciated with field-induced spin disorder and hysteresis in 
these materials. By performing systematic MR measure- 
ments in lightly-doped Ndi.g75Ceo.o25Cu04, which have 
essentially the same AF structure as Nd2Cu04 but with 
enough charge carriers for resistance measurements, we 
can directly compare the transport data with neutron 
scattering results. The following Section will describe 
such comparison. 
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temperatures. The insets show the detailed processes of field- 
temperature hysteresis. 



IV. MAGNETORESISTANCE RESULTS 

In the work of Lavrov and co-workers on 
Pri^gLao.yCeo.oiCuO^S, the similarities between 
the critical field for the spin-flop transition and the rapid 
increase of MR have been taken as evidence for spin- 
charge coupling. Because there exists three spin phases 
in Nd2CuC>4 whereas only one in Pri. 2 9La .7Ceo.oiCu04, 
one would expect some new phenomena related to those 
phases and the transitions between them. Figure 8 shows 
the c-axis MR effect AR c /R c (0) of Ndi. 9 75Ce .025CuO 4 
at different temperatures. As a function of increasing 
field along the [110] direction, R C /R C (Q) initially in- 
creases quickly but then descends slightly. The increase 
in the critical field for spin-flop transition, Bsf, with in- 
creasing temperature for T below 30 K is consistent with 
earlier neutron scattering experiments on Nd2Cu04 14 . 
The increase of f$sF with increasing temperature ends 
above T 3 = 31 K, where the type-Ill to type-II spin- flop 
transition occurs [Fig. 8(a)]. Compared with the results 
on Pri.2gLao.7Ceo.oiCu04^, the data suggest that the 
changes in MR effect originate from the differences of 
TSsf in three different spin phases. Below 30 K, Hsf 



FIG. 8: AR C / R c {0) as a function of increasing magnetic field 
at different temperatures around (a) 30 K and (b) 70K. Note 
that these temperatures are close to Tz and T2, respectively. 
The applied magnetic field is in the Cu02 afc-plane along the 
[1 10] direction. The arrows in (a) indicate the temperature 
dependence of Bsf- The inset in (b) plots the maximum of 
AR C / R c (0) as a function of increasing temperature. It shows 
1/T behavior. Note that the vertical axes are log scale in (a) 
and linear scale in (b). 



increases slightly with increasing temperature. It then 
decreases with increasing temperature beyond 30 K. 
Finally, when the system changes to type-I phase above 
67 K, the sharp increase of R c /R c (0) at low fields almost 
disappears, and all AR C / R c (0) data nearly fall into one 
curve [Fig. 8(b)]. Similar phenomena are also found for 
-R[iio]/-%io](0) and i?[ioo]/#[ioo](0)- The temperature 
dependence of the maximum of AR C / i? c (0) is shown 
in the inset of Fig. 8(b), which can be fit by the 1/T 
function similar to low temperature intensity changes in 
NdaCuO^ and Ndi. 85 Ceo.i 5 Cu04^. 

In previous work, a fourfold angular oscillation in MR 
effect has been identified in Pri.2gLao.7Ceo.oiCu04 19 and 
Pri.ssCeo.isCuO^ for an applied magnetic field in the 
CuC>2 plane. For lightly doped Ndi.g75Ceo.o25Cu04, we 
expect to observe similar fourfold oscillation behavior. 
In addition, we hope to determine whether transitions 
across different spin phases affect the MR. To study the 
angular dependence of resistance, we rotated the sample 
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FIG. 9: Angular dependence of R[uo] at (a) T — 20 K in 
type-Ill phase and (b) T — 35 K in type-II phase, where [110] 
indicates the current direction. The horizontal axes represent 
the angle (in degrees) between the magnetic field and the a- 
or b-axis which cannot be distinguished for the tetragonal sys- 
tem. Note that data in (a) at 1-T data show non- observable 
MR effect while there are clear MR effects at 1-T in (b). 



with respect to the c-axis which remained perpendicu- 
lar to the magnetic field, i.e. the field rotated within 
the a6-plane. The rotation angle is defined to be zero 
when B||[100]/[010]. As shown in Fig. 9(a), a fourfold 
feature in iirno] similar to the earlier results is found as 
a function of the rotation angle. While R\no] along the 
[110]/ [110] direction is much higher than that along the 
[100]/[010] direction at low fields, the situation reverses 
itself at high fields. 

Assuming a strong spin-charge coupling, one can un- 
derstand the microscopic process of the MR effect as 
follows. When an in- plane magnetic field is applied 
along the [110]/[110] direction, spins in the type-Ill 
noncollinear structure of Fig. 1(c) rotate continuously 
to form the collinear structure perpendicular to the 
fieldii While these diagonal (Cu-Cu) directions are easy 
axes in the collinear spin structure with relatively small 
B SF [110], a perfectly aligned field along the [100]/[010] 
direction induces a first-order spin-flop transition with 
a much larger critical field Bs F [100]it For a magnetic 
field in the intermediate directions, it first induces a tran- 



FIG. 10: (a) Angular dependence of R [10 o] at 20 K and 10-T. 
While we still observe the fourfold feature, the MR effects do 
not have the expected 90° symmetry. We speculate that this 
is due to the small misalignment of the sample with respect to 
the applied field, (b) Hysteresis behavior of R[ioo] around the 
transition between the type-II and III phases at 5-T. At each 
measured temperature, a constant has been subtracted from 
the resistance values in processes a and b for clarity. The 
solid and dotted lines represent increasing (b) and decreasing 
(a) temperatures, respectively, as shown in the inset. 



sition into the collinear state and then smoothly rotates 
the spins to positions perpendicular to the fielcUii^. If 
an applied field is less than 1-T, there is no spin-flop tran- 
sition at any field orientation and thus no MR effect, con- 
sistent with Figs. 8(a) and 9(a). When an applied field 
(2-T < B < 6-T) is larger than B SF [Il0] but smaller 
than Bs F [100], one would expect to observe a spin-flop 
transition, and therefore the MR effect, for fields along 
the [110]/[110] direction but not along the [100]/[010] 
direction. This is exactly what we find in Fig. 9(a). Fi- 
nally for B > B sf [100], the aligned collinear spins will 
simply follow the rotation of the field in all directions. 
Here the magnitude of the MR_effect for B||[100]/[010] 
is larger than that for B||[110]/[110], therefore causing a 
new four-fold MR oscillation with 90° angles shift from 
that in B SF [110] < B < B 5F [100]. 

As the temperature of the system is increased to 35 K 
in the type-II phase, the magnitude of Bg F [110] becomes 
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FIG. 11: (a) i?2 — -Ri during increasing and decreasing tem- 
perature processes, as shown by the arrows. The positions 1 
and 2 are defined in the inset of (a). With increasing field, the 
widths of hysteresis become larger. (b)Ri nc — Rdec, the dif- 
ferences between increasing and decreasing temperature pro- 
cesses at positions 1 and 2 are plotted by lines and symbols, 
respectively. 



smaller than that at 20 Ki^i^iS. This explains the ob- 
served fourfold MR oscillations in Fig. 9(b) at 1-T, while 
similar oscillations are only seen for fields above 2-T at 20 
K [Fig. 9(a)]. On increasing the applied field from 1-T, 
the oscillations change from sinusoidal shape to a flatfish 
(or concave) top but the symmetry of the oscillations as a 
function of rotational angle remains even up to the max- 
imum field of 14-T. The data are clearly different from 
that at 20 K for fields above 6-T, where MR effects are 
maximal along the [100]/ [010] directions for B > 6-T. To 
consistently interpret the MR data of Figs. 9(a) and (b), 
we speculate that Bsi?[100] is larger than the maximum 
applied field of 14-T around 30 K. While this scenario 
is difficult to prove because neutron experiments in the 
B||[100]/[010] geometry have not yet been carried out at 
this temperature 17 , a large (> 4.5 meV) in-plane spin- 
wave gap associated with the Nd 3+ -Cu 2+ interactions 
has been reported in the type-II phase of Nd 2 Cu04 3 °. 
If closing such a spin-wave gap is required to induce the 
spin- flop transition in the B||[100]/[010] geometry, the 
critical field necessary to produce Zeeman energy larger 



than 4.5 meV will exceed 14-T assuming only Cu 2+ con- 
tributions [see equation (2) below]. Alternatively, one 
might imagine that the bigger low-temperature MR ef- 
fect in the B 1 1 [100] / [010] direction is somehow related to 
the larger Nd 3+ moments and/or the first order nature 
of the spin-flop transition in this directioniSiSSi. While 
how MR is affected by the Nd 3+ -Cu 2+ coupling is un- 
known, a small misalignment of the sample with respect 
to the magnet around [100]/[010] directions can affect 
dramatically the observed MR. Such misalignment may 
also explain the slightly different MR values at 45° and 
135° in Fig. 9(b). We note that similar, but less obvious, 
behavior is also present in Fig. 9(a) and in previously re- 
ported MR datai2*2fi. 

Figure 10(a) shows the angular dependence of i?[i o] 
at 20 K and 10-T. While i?[ioo] has fourfold oscillations 
similar to that of R[no] at the same temperature [Fig. 
9(a)], the MR differences between 45° and 135° are more 
obvious. With increasing temperature, the fourfold os- 
cillations are replaced by a square-wave-like feature [Fig. 
10(b)] similar to Fig. 9(b). Note that our neutron scat- 
tering revealed a clear hysteresis through the type-Ill to 
II collinear phase transition [Fig. 7(a)]. To see if MR fol- 
low such hysteresis, we performed careful measurements 
on field-warming and cooling as processes b and a, re- 
spectively [see the inset of Fig. 10(b)]. The outcome in 
Fig. 10(b) shows clear hysteresis across the type-Ill to II 
transition, consistent with the neutron scattering results 
of Fig. 7. In addition, we find that the relative value of 
i?r 100 i shifts 90° across the transition, suggesting that the 
MR effects are sensitive to the differences in the type-Ill 
and II collinear spin structures [Fig. 1(c)]. If we define 
the resistance at 135° position as i?2 and 225° as R2 
is larger than R\ in the type-Ill collinear phase while the 
reverse is true in the type-II collinear phase. 

Using the resistance difference between positions 1 and 
2, we probe the phase transition between type-Ill and 
II collinear states in great detail. If there is no field- 
induced hysteresis, R2 — Ri should be the same for either 
warming or cooling. Figure 11(a) indicates that this is 
not the case. On warming, the type-Ill collinear to II 
collinear transition temperature increases with increas- 
ing field. On the contrary, the type-II to III transition 
temperature decreases with increasing field on cooling. 
As a consequence, the width of the hysteresis increases 
with increasing field and can be as large as ^15 K at 
14-T [Fig. 11(a)]. At 5-T, the width of the hysteresis is 
about 2 to 3 K, completely consistent with the neutron 
scattering results of Fig. 7. Figure 11(b) shows the differ- 
ences in resistance at positions 1 and 2 between increasing 
and decreasing temperature processes for various applied 
fields. The results also suggest an increasing hysteresis 
in the phase transition with increasing magnetic fields, 
consistent with Fig. 11(a). 

Finally, we describe the transport measurements asso- 
ciated with the c-axis disorder seen by neutron scattering. 
Following the same processes as the inset in Fig. 7(a), 
we have measured the resistance of Ndi.g75Ceo.o25Cu04 
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5-T through processes a, b, and d similar to that defined in 
Fig. 7(a), except in this case the final temperature was fixed 
at 34 K. Process a does not across the phase transition tem- 
perature, and therefore does not exhibit disorder. Processes 
b and d have spin disorder. While clear differences are seen 
in R c , there are no observable differences in i?[ 100 i for these 
different processes. 



in three different current directions, R[iqo], -Rfiiob and 
R c . For the in-plane resistances, we find no distinguish- 
able difference between i?[ioo] and R[no] after waiting one 
hour for each measurement. The results of i?[ioo] at 34 
K and 5-T in Fig. 12 show overlapping curves for pro- 
cesses of a, b, and d. On the other hand, R c displays 
clear distinctions among the varying processes as shown 
in Fig. 12. Since the resulting differences in i?[ioo] for 
processes a, b, and d are less than 5 x 10 -5 , we can safely 
conclude that the observed deviations in R c among these 
processes are intrinsic and may originate from the spin 
disorder along the c-axis. More work is needed to un- 
derstand precise relationship between spin disorder and 
charge transport properties. 



V. DISCUSSION 

Until now, the most successful theory to understand 
the spin properties of Nd2Cu04 is based on pseudo- 
dipolar interaction (PDI) originally proposed by Van 
Vleck in 1937, 



(5) 



where I and I denote the lattice sites and the function 
V(R) decreases faster than R~ 3 as R — > oo. To explain 
the reorientation of the spin structure, Sachidanandam et 



alH considered three major interplane interactions be- 
tween Nd-Nd, Nd-Cu and Cu-Cu, labeled as X, Y and Z 
respectively in Fig. 1(a). The interactions X and Z tend 
to generate the type-Ill or I spin structures, and while 
Y prefers type-II phase. The interactions between spins 
are proportional to their local susceptibilities (m) . Since 
mjvd is proportional to 1/T— and mcu varies little below 
40 Kii, X cx 1/T 2 , Y oc 1/T and Z » Constant. With 
decreasing temperature, Cu-Cu (Z) interactions initially 
turn on below T\ and Nd 2 Cu04 orders antiferromagneti- 
cally with the type-I spin structure [Fig. 1(a)]. On cool- 
ing to intermediate temperature T%, Nd-Cu (Y) inter- 
actions become important and the system transforms to 
the type-II noncollinear spin structure. Finally below T3, 
Nd-Nd (X) interactions dominate and induce the type- 
Ill noncollinear spin structure (Fig. 1). 

The noncollinear spin structures of Nd2Cu04 have a 
small spin-wave anisotropy gap Ao at zero fieldM. When 
an in-plane field is applied, the Zeeman energy shifts 
the spin-wave dispersion and closes the anisotropy energy 
gap, resulting a transition from noncollinear to collinear 
spin-flop phase. Petitgrand et alJ& have given the crit- 
ical field of spin-flop transition when the field is along 
[1101 direction, 



B SF [no] = 



A 
gma B 



(6) 



where Ao is in-plane spin-wave gap at zero field, g is 
Landau factor, m the effective moment, and ug the 
Bohr magneton. This equation has been successfully 
used to explain the temperature dependence of Bg^ [110] 
for Pr 2 Cu04i&. For Nd 2 Cu0 4 , the in- plane Cu spin- 
wave gap has a 1/T dependence and the out-of-plane 
gap is essentially temperature independent^. In ad- 
dition, Nd spin-waves exhibit anisotropic gaps at low 
temperatures^ 2 -. Since an applied field of a few Tesla in 
the Cu0 2 will not change the large (> 5 meV) Cu spin- 
wave gap in type-Ill phase below 30 K—, the spin-flop 
transition there is most likely induced by closing the Nd 
spin- wave gap. 

The existence of a spin-charge coupling has been sug- 
gested in lightly electron-doped Pr1.29Lao.7Ceo. 01 C11O4, 
but a detailed microscopic understanding of how such 
coupling occurs is still lacking 1 ^,. If the itinerant 
electrons are coupled to the localized spins directly, 
one would expect to observe their signatures in the 
zero-field resistance when Pri. 2 gLao.7Ceo.oiCu04 and 
P r i.85Ceo.i5Cu04 order antiferromagnetically 19,20 , and 
when different noncollinear spin phase transitions oc- 
cur in Ndi.975Ceo.o25Cu04 (Fig. 2). However, AF or- 
der appears to have no observable effects on zero-field 
resistance. Instead, the maximum of AR C / i? c (0) in 
Ndi.975Ceo.o25Cu04 shows a 1/T temperature depen- 
dence [inset of Fig. 8(b)], very similar to the 1/T 
temperature dependence of the Nd moment in various 
Nd-containing Nd 2 _ x Ce ;c Cu04 compounda 11 ! 14 ! 15 ! 18 ' 29 . 
This strongly suggests that the observed MR effects in 
Ndi.975Ceo.o25Cu04 and other electron-doped materials 
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are somehow related to the rare earth (Nd,Pr) moments 
and/or Nd(Pr)-Cu coupling. This picture may also ex- 
plain why, when the dominant spin-spin interactions are 
from Cu-Cu with negligible Nd moments in the type-I 
collincar phase (T > 68 K), the weak MR data are essen- 
tially temperature independent and collapse onto a single 
curve [Fig. 8(b)]. 

To compare our results with hole-doped materials, 
we note that large anisotropic MR effects have already 
been reported for lightly doped La2_^Sr x Cu0^i& and 
YBa 2 Cu 3 Oe+ x 33 i 34 . These results have been interpreted 
as due to the influence of an applied magnetic field on 
stripesiSiSi 3 ^ S pi n _ rbital coupling 35 , redistribution of 
magneto-elastic antifcrromagnetic domains 3 ^, or canted 
AF spin structures^. At present, there is no consensus 
on a microscopic picture for the MR effects in hole-doped 
copper oxides and more work is needed to test the pre- 
dictions of different models. However, regardless of the 
details for each model, what is clear is that transport 
properties of electron or hole doped copper oxides are 
closed related to the AF order in these materials. 



electron-doped copper oxide affects both the in-plane 
and out-of-plane MR. The application of a in-plane mag- 
netic field can induce c-axis spin disorder and hystere- 
sis in the AF phase transitions. By comparing neutron 
scattering results of Nd2CuC>4 with the MR effects in 
Ndi.975Ceo.o25Cu04, we show that the transport prop- 
erties of these materials are very sensitive to the subtle 
changes in the spin structures. Our results thus provide 
further evidence for the existence of a strong spin-charge 
coupling in both electron and hole doped copper oxides. 
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